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Succeeding in predicting + and 1 + states of D and D s heavy mesons with our semi- 
relativistic quark potential model, we examine a method for constructing Lorentz-invariant 
scattering amplitudes and/or decay widths and develop a formulation to calculate Lorentz- 
qq ' boosted ones given the rest-frame wave functions in our model. 

To illustrate the effectiveness of our method, we apply the formulation to the problem of 
^/-^ . calculating the semileptonic weak form factors from the rest-frame wave functions of heavy 

mesons and numerically calculate the dynamical 1/mQ corrections to those for the process 
00 ! B -> D M £p based on our model of heavy mesons. It is shown that nonvanishing expressions 

. for pi(u>) = p2{uj) and Ps{uj) = pi{uj) = are obtained in a special Lorentz frame, where 

Pi(u>) are the parameters used in the Neubert-Rieckert decomposition of form factors. Various 
' values of form factors are estimated. These values are compatible with recent experimental 

, data as well as other theoretical calculations. 

o 

-p^ ' §1. Introduction 

The discovery of the narrow meson states D s0 (2317) by BaBar 11 and £> s i(2460) 
by CLEO 2 ' and the following confirmation by Belle® ® has motivated many theo- 
^ . rists to explain these states, as previous studies of these states (Refs. [B]) and [7j)) 

using quark potential model apparently fail to reproduce these mass values. More 
recent experiments have found many other heavy-light mesons : broad Dg(2308) 
and 1)^(2427) mesons found by the Belle collaboration,® which are identified as 
cq (q = u/d) excited (I = 1) bound states and have the same quantum numbers, 
j p = + and 1 + , as D s j, respectively (The masses of these mesons have been re- 
ported with quite different values in the CLEO and FOCUS experiments, and hence 
are not yet determined with definitenessP'^); narrow B and B s states of I = 1, 
Bi (5720), 5|(5745), and £* 2 (5839) found by CDF and DO, 9 ' whose decay widths 
are also narrow because of their decay through the D-waves; and seemingly ra- 
dial excitations (n = 2) of the + -D<,o(2860) state found by BaBar-^ and the 1~ 
Ds* (2715) state found by BelleP^ Furthermore several cc quarkonium-like states 
have been discovered : X(3872), X(3940), y(3940), Z(3930), and y(4260). 

The discovery of these mesons triggered a series of studies of the spectroscopy 
of heavy-light and/or heavy-heavy hadrons using many kinds of ideas, including our 
semirelativistic quark patential model. Because we do not discuss these ideas in this 
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paper, we refer the reader to the review articles .OJEDJID i n previous papers,® we 
formulated a method for calculating the spectrum of heavy-light mesons in order to 
construct the Schrodinger equation for a bound state. In that method, the mass is 
expressed as an eigenvalue of the Hamiltonian of the heavy-light system, in which 
a bound state consists of a heavy quark and a light antiquark, and negative energy 
states of a heavy quark appear in the intermediate states when calculating an energy 
eigenvalue. In oder to make our formulation more reliable, we attempt to reproduce 
narrow B and B s states of £ = 1, Si (5720), B|(5745) and S* 2 (5839), together with 
higher states of J p for D, D s , B and i? s ,® an d also to reproduce radial excitations 
of the 0+ state of £> s0 (2860) and the 1" state of Ds* (2715). ^ These calculations 
show our model's reliability, as these results fit well with experiments. What we need 
to do next is to show that our approach can also produce a method for calculating 
scattering amplitudes and decay widths using the rest-frame wave functions. Here, 
we concentrate on the construction of such a formulation, and therefore we refer 
readers the review paper Ref. I18|) for explanations of other methods and their results, 
because there are many papers that give estimates of widths of many decay modes 
of heavy mesons. 

The heavy quark effective theory (HQET) is a very convenient tool for studying 
heavy mesons/baryons, including at least one heavy quark.® >EB T n particular, there 
are many studies on the Isgur-Wise function ED~E3 j-,^ m ost Q f the calculations in 
those works use a simple-minded Gaussian form for the meson wave function or a so- 
lution to a single-particle Dirac equation in a potential. Some other people have used 
the Bethe-Salpeter solution for a meson wave function, but it seems there has been 
no systematically developed way to calculate the higher-order terms in 1/toq®'® 
although some papers® claim to have done so. Other people® used the light-front 
formalism to calculate the scattering amplitudes, so that in the heavy quark limit 
they obtained the Isgur-Wise function, in which they adopted the Gaussian-type 
wave function for the heavy mesons. There is a recent paper® in which the authors 
use a " relativistic" formulation to calculate the masses of heavy mesons and apply 
it to calculate the Isgur-Wise functions for semileptonic B decays. Although that 
paper adopts an approach that is different from ours, they present a method for 
calculating higher-order corrections in I/tuq to form factors, which is one of the 
purposes of this paper. 

Most of the other calculations of form factors used the so-called trace formalism 
developed by Falk et al.p® which assumes the Lorentz covariance of the final form 
and fits it to the experimental data, or the QCD sum rule® was used to obtain 
the relation with the quark/gluon condensates. In previous papers,® we developed 
a semi-relativistic formulation for calculating the meson masses and wave functions 
for the heavy- light Qq system with one heavy quark, Q, and one light antiquark, q, 
introducing phenomenological dynamics. That is, assuming a Coulomb-like vector 
potential as well as a confining scalar potential to Qq bound states, we expanded the 
effective Hamiltonian and then perturbatively solved the Schrodinger equation as an 
expansion in 1/mQ. The meson wave functions thereby obtained and expanded in 
1/mg can be used in principle to calculate ordinary form factors or Isgur-Wise func- 
tions and their corrections higher-order in 1/ txiq for semileptonic weak or other decay 
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processes. However, what we have obtained are wave functions in the rest-frame, 
and therefore we need to develop a method to obtain Lorentz-invariant amplitudes or 
Lorentz-boosted wave functions. In this paper we would like to address the problem 
of how accurately we can calculate scattering amplitudes and/or decay widths of 
heavy-light systems using the rest-frame wave functions. In the following, explain- 
ing the results obtained in Ref. [T5j) . we show how to calculate the Lorentz-invariant 
amplitudes by taking the Isgur-Wise function as an example. This calculation is 
general enough to apply our formulation to other scattering/decay amplitudes. 

Following a previous study, 35 ' in which we derived the zeroth-order form of the 
Isgur-Wise function expanded in 1/toq, we develop a formulation for calculating 
higher-order corrections to the semi-leptonic weak form factors in 1/mg from our 
semi-relativistic wave functions for heavy mesons in the rest-frame obtained in Ref. 
[T5|) . The problem involved in constructing a Lorentz-boosted wave function is that 
there is ambiguity in determining the space-time coordinates of two composite par- 
ticles from information regarding one bound state, in this case a heavy meson. We 
study four cases for the reference frames of composite particles and then give a 
prescription for calculating matrix elements of currents using the rest-frame wave 
functions in In £j3j we give the form factors at zeroth order for four reference 
frames, show that they agree with each other in the HQET limit, and then obtain 
them at first-order for one special reference frame. The calculation to lowest order 
in 1/rriQ gives the numerical value of the slope for the Isgur-Wise function at the 
origin, and the semi-leptonic weak form factors are calculated up to first-order in 
1/mQ in It is found that there are no dynamical contributions to the form fac- 
tors; i.e., the first-order corrections to the wave functions do not contribute to the 
form factors. In SjH studying B — * Dtv and B — * D*lv processes, physical quantities 
related to the CKM matrix elements are obtained. Finally, £}5]is devoted to summary 
and discussions of the results obtained in this paper. All the details associated with 
calculations of the form factors are given in the appendices. 

§2. Formulation 

2.1. Schrodinger equation 

In Ref. [T5|) . we calculated the mass spectrum of heavy mesons in a so-called 
Cornell potential model, which includes both scalar and vector potentials. Although 
we needed only a stationary system, or rest-frame wave function there, we must now 
treat a moving frame of a heavy meson to calculate weak form factors. Although 
we cannot work in a fully relativistic system as far as the Hamiltonian formulation 
is concerned, we adopt the notation of the Nambu-Bethe-Salpeter equation as much 
as possible. First, we define a heavy meson (A) wave function composed of a light 
anti-quark, q c (x), and a heavy quark, Q(y), as 

<0| q c a (t, x) Q p (t, y) | A; P x {£}) = tf Xafi ((0, x - y) ; P x ) e" lP * * (1) 

where a and (5 are Dirac indices of four-spinors, Px^ is the heavy meson four- 
momentum, {£} is a set of quantum numbers, and Px ■ y = Px^y^- Note that this 
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definition of the wave function, i/j x , equals to that used in Ref. [T5]) multiplied by -1. 
The heavy meson state is an eigenstate of the four-momentum operator given by 



V f ,\X;P x {£})=P x ^\X;P x {£}), P Xo = ^M x + P x , (2) 

with M x a heavy meson mass. In this case, the Schrodinger equation in the X 
moving frame is given by 



Hip x ((0, x-y)-P x ) = + P| ((0, x-y);P x ), (3) 

where the Hamiltonian is given by 

H = (a q ■ p + (3 q m q ) + a Q ■ (P x - p) + P Q m Q ^ +^2(3 q O q iVij(x - y)(3 Q Qj . (4) 

Here, p = —i'Vx and P x is the heavy meson momentum. In our case, Eq. (|4|) gives 
O q i = Oqi = 1 for the scalar potential, Vij{x — y) = S(r) and O q ; L = ^ q/i and 
OQj = ^q 11 for the vector potential, Vij(x — y) = V(r) with r = \x — y\. Actually, 
we have included the transverse part of the vector potential and the total potential 
in Eq. ([!]) is given by 

f3 q (3 Q S{r) + jl - ^ [d q ■ d Q + {a q ■ n) {a Q ■ n)]| V{r), (5) 

with 

T A- Oil V 

S(r) = - + b, V(r) = --^, n=~. (6) 
a z 3 r r 



In Ref. [15]) . we solved the Schrodinger equation given by Eq. (J3|) with P x = order 
by order in 1/mg by expanding H, tp x , and M x in 1/mg, in order to numerically 
calculate the mass spectrum, i.e., an eigenvalue, M x , in each order of perturbation, 
and then to relate the results to those of HQET. 

2.2. Lorentz boost 

In order to calculate the decay rate, we need to Lorentz boost the rest-frame 
wave function that we obtained in previous papers. To explain the idea, we consider 
the B — > D^lv decay process as an example and express the wave function in the 
moving frame in terms of that in the rest-frame. 

Assume that the heavy meson is moving in the +z direction with a velocity 
/?. Below, we consider two cases of composite particles, i.e., q c and Q. In the 
following, for the initial state of a heavy meson, we adopt the notation q c (x°, x) and 
Q(y°,y) and for the final state, q c (x'°,x') and Q(y'°,y'). The relation between the 
constituent particles in the rest and moving frames is given by 

Gq^x^^G-^G-lq^x' ,?'), (7) 
QQ^y^^Q- 1 = G~lQ p {y'\y'), (8) 
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where Q is the Lorentz-boost operator, G is its spinor representation, and their 
definitions are given by 

X;(m x ,0){£}) = \X;P x {£}) (9) 



G | 

G = cosh — + q 3 sinh — , /3 = tanhw, 7= — = cosh (p. (10) 

2 2 ^> r ^ ^ > 

We have assumed that the heavy meson is boosted in the +z direction. Here, we 
have a = 7°7, and f3 and 7 are the velocity and Lorentz factor of the heavy meson, 
respectively, which can be expressed in terms of one parameter, ip. 

Even if two constituent quarks have equal times in one frame, they have different 
times in another frame. In order to pull back the time of the heavy quark to that of 
the light quark, we make an approximation that the heavy quark propagates freely 
over a short time interval, assuming 

Qp {t + 8t, v) - ex P (-imQ-ySt) Qp {t, y) . (11) 

Below we consider two cases. In the first case, the two constituent quarks have 
equal time in the rest-frame of the meson, and we apply the above approximation 
to the heavy quark in the moving frame. In the second case, the two constituents 
have equal times in the moving frame of the meson, and we apply the approximation 
to the heavy quark in the rest-frame. In these cases, different relations are derived 
between the wave function of the meson with finite momentum and that of the meson 
at rest. These are as follows : 

i) t = x° = y° and x'° 7^ y'° (equal time in the rest-frame) 

x = x cosh ip + t sinh ip, x = t cosh cp + x sinh (p, 

y =y cosh ip + 1 sinh ip, y =tcos\np + y sinh(/?. (12) 

In this case, an approximate relation between the rest and moving frame wave func- 
tions is given by 

^ a/3 ((0,x);Px)^G Q7 G^^ 7 ,((0,f ± ,7- 1 x 3 );(M x ,o)) e< M *~ m 

(13) 

ii) t' = x'° = y'° and x° 7^ y° (equal time in the moving frame) 

x = x cosh ip — t sinh (p, x = t cosh ip — x sinh ip, 

y 3 = y' 3 cosh ip — t' sinh ip, y° = t' cosh p — y' 3 sinh (p. (14) 

In this case, an approximate relation between the rest and moving frame wave func- 
tions is given by 

^ Q/3 ((0,i ? ); J Px)^G Q7 G /35 Vl 75 ((0,x ± ,7x 3 );(M x ,0)) e K M x- mQ ),^\ (15) 



The derivation of the Lorentz boosted wave functions given by Eqs. (fTBj) and (|T5|) is 
given in Appendix [Bj 
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2.3. Wave function 

The explicit form of the wave function we obtained in the previous paper has 
the form 

Ux{p)^x (x) = i> l x FWT (x), (16a) 
U x (p) = U c U FWT (p), v4(^) = V>x((0,x);(m x ,o)), (16b) 

^xfwt(^) = 14c o(x) +tfx\(x) + ■■■. ( 16c ) 
¥xo@) = = (o #} m ) ■ (16d) 

We have solved the Schrodinger equation in terms of the wave function ip x fwt > n °t 
if) x . Here, the 4x4 wave function, ip x , is transformed just once with the Foldy- 
Wouthuysen-Tani transformation, C/fwT: acting on the heavy quark and with the 
charge conjugation operator, U c , into V'xFWT' an d ^ s *~th ° r der term in 1/mg is 
given by ip Xi - Higher-order corrections depend on the positive and negative compo- 
nents of the heavy quark, Sl/f 1 , which are given by 

W+W = VMx~ (o ^ m (x)) , Vfffl = TMx (^ m (x) 0) . (17) 

Here, p = pq appearing in the argument of Ux(p) above is the initial momentum of 
the heavy quark, and henceforth, the color index (iV c = 3) is omitted, because the 
form of the wave function is the same for all colors, £ stands for the set of quantum 
numbers, j, m, and k, and the positive/negative component wave functions, , are 
given in terms of a 4 x 2 wave function, \Pj m {x), as 

■?-<*> = ?( -i^rUft.g) ) < 18 > 

where r = \x\, n = x/r, j is the total angular momentum of the meson, m is its 
z component, k is a quantum number that takes only the values k = ±j, ±(j + 
1) and ^ 0. (More details are given in Appendix ICl) 

We showed in the previous paper that &£{x) is an eigenstate of the operator 
K = —f3 q {S q ■ i + 1), with eigenvalue k, and we classified the spectra in terms of 
k. People normally classify the spectra in terms of sf , which is defined to be an 
eigenvalue si [sg + 1) of the operator (s^) 2 ,® together with parity, P. We can show 
the following direct relation among k, se, and P; i.e., sg and P can be explicitly 
written only in terms of k as^ 

k = ± (s t + ^ or s e = \k\ - i (19a) 

P = ^(-l) |fc|+1 , (19b) 
where we have definitions 
s e = 1+ ±E q , (s e ) 2 = V + £■ U q + |, /= -if x V, 4 = ( o ^ ) • ( 2 °) 
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The corresponding physical quantities between k and S£ and the first few states are 
listed in Tables [J and [TH respectively. 

Table I. Corresponding values for k, si, and P. 



p 




(-Y 


(-Y 




k 


-(i + i) 


3 + 1 


-3 


j 


St 


i + i 


3 + i 


i-i 






-i 






-1 



Table II. States classified with respect to various quantum numbers. 



■1> 

3 


0~ 


1" 


0+ 


1 + 


1 + 


2+ 


1" 


2~ 


2" 


2+ 


k 


-1 


-1 


1 


1 


-2 


-2 


2 


2 


-3 


3 


sf 


1 - 

2 


1 - 

2 


1 + 

2 


1 + 
2 


3 + 

2 


3 + 

2 


3 - 
2 


3- 

2 


5 - 
2 


5 + 
2 






K l 










3? 









From Tables U and HH it is seen that there are both advantages and disadvantages 
of using these quantities to classify the heavy mesons. The quantity sf has an 
intuitive meaning that it represents the total angular momentum of the light degrees 
of freedom, while k does not. However, k includes information concerning parity, 
while si does not. Hence we need to add P, writing it as sf. Certainly either one can 
be used to classify the heavy meson spectra equally well; i.e., they are equivalent. 

2.4. Normalization 

Normalization of the state \X;Px {£}} is formally given by 

(X; P' x {£'} \X; P x {£}} = (2tt) 3 2P Xo 5 3 (P x - P' x ) 5 e ,t (21) 

or 



afz tr 



tf? ((0, z) ; Px) A ((0, z) ■ P X ) ~ 2P X0 b it 



(22) 



whose derivation is given in Appendix [Bl Actually, the normalization given by 
Eq. (|22p does not hold, because of approximations adopted in the previous subsec- 
tion, and hence, instead of Eq. (|22p . we define the normalization in the rest-frame 
as follows : 



d 3 z tr 



Mx,0 h& (0,2); M X ,0 



2M x S t 



(23) 



In order to calculate the normalization in a moving frame, we have to specify the 
quantum numbers in which we are interested. Because we would like to compute 
the Isgur-Wise function and its higher orders, the cases of the pseudoscalar state 0~ 
and vector state 1~ are calculated in the following. 

The normalization condition given by Eq. (|23p can be rewritten in terms of 

V'x fwt( j; ) as 



d 3 xtr 



^xfwt* (z) ^xfwt (£) = 2M x Sa'. 



(24) 
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All the details concerning the wave functions are given in Appendix ICl 

In order to calculate the normalization of the wave function in a moving frame 
and/or the matrix elements of the semi-leptonic decay processes, we need to develop 
a formulation of their calculations in terms of the rest-frame wave functions. Let 
us assume that a physical quantity is already given in terms of the rest-frame wave 
functions and rewrite them in terms of the FWT-transformed ones as 

d 3 x (O q ) aa , (O q )^, A' a %B'^? jS A a/Y Bp S ,^ Xj , s , 
= J d 3 xtr (y t xU'iO q Atf x B' F OlB'*) 

= j d 3 xtr [^ FWT A'^O q Ai; XFWT U x 1 (PQ)B T OlB'*U x/ (p' Q )] , (25) 

where (O q ) aa i, A' ai , and A a y act on light quarks, while (Oq B'^g and Bps' 
act on heavy quarks, and use has been made of 

= U x x (jpq) <8> iPxfwt(pq) = V'xfwt V^ipo), i> x , = 4'fwt u x 1 (p'q)i 
U x \p Q )B^OlB'*U xl {p' Q ) = U XFWT ( PQ )U- l B T OlB'*U c U x , FWT (p' Q ). (26) 

In the course of obtaining relativistic results from the rest-frame wave functions, 
we calculate the normalization of the wave function in a moving frame in two extreme 
cases, i) t = x° = y° and ii) t' = x'° = y'°. The FWT-transformed rest-frame wave 
function is given by, to first order in l/mg,® 

^fwt(O-) = *±i + c^' + O (l/m|) , (27a) 
^fwt(11 = + cr+' 2 ^2 + + <h- + Ci-~ 2 *-2 + O {l/m 2 Q ) , (27b) 

where only the value of the k quantum number is written as a subscript of 
The total angular momentum, j, though not given explicitly, should be the same 
on both sides of the equations and the coefficients, c 1 ! t , are given in [T5j) . Their 
explicit expressions are not necessary here, since we show that the higher-order 
corrections do not contribute to the physical quantities at this order, 0(1/ mq). 
Actually, Eq. (|27p can be derived from the conservation of total angular momentum 
and parity alone, without explicit interaction terms specified, since we know the 
complete set of eigenfunctions with j, m, and k quantum numbers. (Details are 
given in Appendix O) 

We calculate the normalizatio up to first-order in 1/mn below. 

i) t = x° = y° and x'° ^ y'° (equal time in the rest-frame) 

Using Eqs. (fT3j) . (f23| and (f25|) . we obtain, to first-order in 1/mg and both for 0~ 
and 1~ states given by Eq. ([27|) . 

d 3 x tr ((0, x) ; P X ) tf x ((0, x) ; P x )] = 2M x7 3 + O (l/m 2 Q ) . (28) 

ii) t' = x'° = y'° and x° ^ y° (equal time in the moving frame) 
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Similarly using Eq. (|15p . (I24j) and (|25h . we obtain 



V£ + ((0, x) ; Px) tf x ((0, a?) ; P x )} = 2M xl + O {l/m%) 



(29) 



cTa; tr 

Only this case agrees with the relativistic normalization given by Eq. (|22p . 

2.5. Matrix elements 

The light anti-quark is treated as a spectator; that is, it does not interact with 
other particles, except for gluons represented by potentials in our model. The heavy 
quark has a weak vertex, and its current is in general given by 



j(t,x) = Q x ,Jt,x) O a/ 3Q XfS (t,x), 



(30) 



where O is a 4 x 4 matrix. By inserting the number operator of the anti-quark and 
by assuming vacuum dominance among intermediate states, we obtain the formula 
for the matrix element of the above current, Eq. (|30p . between heavy mesons : 



(X';P X , {£'}\j(t,x)\X;P x {£}) 

l 3 y tr [^ x , \(0,y-x)-,P x ,)(o®4> x ((0,y-x);P x : 

-i{P x -P x ,)-x 



x e 



(31) 



or 



X';P X >{£?} j(0,0) X;P X {1} 



~ / d x tr 



^((0,5)^x0 i> x ((0,x);P x ) O 1 



(32) 



Now we would like to evaluate the matrix element of the process m 
two frames, the B rest-frame and the Breit frame. (The velocities of B and D have 
the same magnitude and have the opposite directions.) In each frame, we apply two 
approximate relations, Eqs. (|12|) and (|14|) . As mentioned in the previous subsection, 
the normalization condition Eq. (|22p of the wave function holds in the approximation 
ii) (equal time in the moving frame) but not in the approximation i) (equal time in 
the rest-frame) . In the approximation i) , we need to renormalize the matrix elements 
by normalizing the wave function, and we redefine the matrix elements of the current 
as follows : 



£>(*) 



j(o,o) 



B 



^2P Bo 2P dMq J d 3 x tr U e D l ((0, x) ■ P DW ) ^((O, x) ; P B )0^ 



(33) 



fd*x tr rP £ '((0,x);P DM ) J d?x tr ^ § ((0,x) ; P s ) 



The numerator here is calculated below using the results of ^2.21 while the denomi- 
nator is calculated in §2.41 and given by Eqs. (|28|) and (|29p . 
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1) B rest- frame 

1-i) t = x° = y° for the £>W meson 

In this case, the matrix element is given by, with our approximations, 



/ 



d 3 x tr 



= J ci 3 xG; 5 G* e ^ w * e ((0,x ± ,7-V);(M DW ,0))o Q ^| 7 ^(^) 

x e ~ i ( M D(*)~ mc )' r l 3x3 

1-ii) t' = x'° = y'° for the meson 

In this case, the matrix element is given by 



d 3 x tr ((0,x);P DW ) ^((0, x) ; Pg)C J 

= J d 3 xG; 5 G; e ^ w ; e ((0,x±,7x 3 );(M DW ,0))o a/3 4 7/3 ( 3 

x e -<^ D (»)-^c)7/3^ 3 _ 

2) Breit frame 

2-i) t = x° = y° for the B and mesons 

In this case, the matrix element is given by, with our approximations, 



d 3 x tr 



((0,x);P DM )^ s ((0,x);P B )O i 



I 

= 7 J <Pxifa.)' 1 JQ(&0(T 1 ) a ^ p W 



x cxp 



- - (M DM — m c + Mb — m b ) 7 2 (5 x 3 



2-ii) t' = x'° = y'° for the B and mesons 
In this case, the matrix element is given by 



d 3 x tr 



^ ( l ) ((0,x);P DW )VM(0,f);P s )O J 



/ 

= 7 _1 / d 3 x^; ia {x)(G^OG^) a ^%^{x) 



x exp 



(M DW -m c + M B - m b ) (3 x 3 



§3. Semi-Leptonic weak form factor 



(34) 



(35) 



(36) 



(37) 



Using the results derived in the previous section, we can now calculate semi- 
leptonic weak form factors including the Isgur-Wise function. In this case the cur- 
rents are given by 

j / , = C t(t,x) 7 ° 7M 6(t,x), (38a) 
jJ = ct(t,£)7Vy 5 6(t,£). (38b) 



Construction of Lorentz- Invariant Amplitudes in HQET 
Define the six independent semi-leptonic weak form factors as^ 1 



11 



(D\j„(p,0)\B) = VM b M d (£+(u)(v b + vd)^ + £-(u)(v B - v D )») , (39a) 



(D*\j^(0,6) \B) = i yjM B M D *Zv(u) 



'pvpa 



* v p 

e Dn, V 



D* U B' 



(39b) 



with 



(D*|jJ(0,6) \B) = ^M B M D * ((l+u)UM < 

-U 2 {u) (e* • v B ) v B p~ U 3 {u) (e* • v B ) v D * ^ , (39c) 



P x /M x , for X = B, 2)W, 



J X ~ * X/ ll ^> iw,^-^, ^ , lo = v b -vd- (40) 

3.1. Zeroth order (Isgur-Wise function) 

After some calculations, all the above form factors are found to be proportional 
to the Isgur-Wise function, or vanish and are given by 

= £ v (u) = Ui(") = U.M = (41a) 
e_(w) = UM = 0. (41b) 

Below we only briefly show how to calculate = from Eq. (j39a|> . since the 

other functions are obtained similarly. All the details are given in Appendix [El We 
define the Isgur-Wise function as 



D 



jo (0,0) 



B 



y/2P BQ 2P D0 J d 3 x tr Vdo ((0, x);P D )^ BQ ((0, a>) ; P B ) 



fd 3 xti\^ DQ ((0,x);P D )\ z fd 3 xtr ^ Bq ((0,x) ;P B ) 
^M B M D i{uj){v B + v D ) . 



(42) 



The expression Eq. (|42j) can be further reduced in two cases of the B frame as follows. 
1) B rest-frame 

In this frame, B is at rest, D is moving in the +z direction with velocity (3 and £(u;) 
is given by 

f d 3 x tr [v£ T o ((0, x) ; P D ) ^ Q (x) 
CM = 



(43) 



y/M^(l+u)JJcPx tr ((0, x) ; P ) 



with the relations 



Pi 



B0 



M B , P do = M d v d = M d uj, 



vdo = v B ■ v D = 7 



1 



y B0 = ^ 



(44) 



2) Breit frame 

In this frame, B is moving in the +z direction, D in the — z with the same velocity 
P/2, and £(u) is given by 

/ d 3 x tr U£ ((°» *0 ; p o) v£o ((°> f ) ; 
eM = 1 J ==, (45) 



/ d 3 xtr|^ ((0,f);P D )| 2 / d 3 ztr ^((0,3);^ 
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with the relations 

P Sq = Mb1, Pdo = M d -/, v So = v DO = -f 



V1-G9/2) 2 ' 



«fi-« D =7 a (i + ^ r )=«, ^ = VhF' f = /^TT- (46) 

We note that here and below we have used the same notation for the velocity and 
Lorentz factor both in the B rest-frame and the Breit frame, although they are 
different quantities, as one can see from Eqs. (JMJ) and (|4Uj) . 

We calculate this Isgur-Wise function in four cases below. The details are given 
in Appendix lE.il 

1-i) B rest-frame and t = x° = y°. In this case, we have 

i{u:) = l-\^El^) + ^+ 0(/? 4 ), (47) 

where 

(r 2 } = J d 3 xti ^xo\x)r 2 ^ X0 (x) , for X = D, or B, (48) 

both of which give the same result, since this is the lowest-order calculation. 

1- ii) B rest-frame and t' = x'° = y'° : 

£M = 1 - \^& D <r 2 } -£+ 0(/3 4 ). (49) 

2- i) Breit frame and t = x° = y° : 

£M = 7' 2 - I (f) 2 l\E B + E D f <r 2 > + 0(/3 4 ). (50) 
2-ii) Breit frame and t 1 = x'° = y'° : 

CM = 7" 2 - I (§) 2 7" 2 (^ + ^) 2 <r 2 > + 0^). (51) 
In the above we have used 

E D = M D - m c , Eg = Mb — m b . (52) 
The rhs of Eq. (|52p can be expanded in 1/rriQ as 

* x = Mx - mQ = ('"« + g (S?) - m ° = g w" (53) 

where the f-th order in 1/mg of Mx is given by C x /(mQ) 1 , and it is shown in Ref. 
[T5]) that depends only on the light quark mass m q , and hence C l D = C B = C l {m q ) 
in our case, since light quarks are either u or d, and we have set m u = m^. Equation 
(|53p is for the 0~ state. When it is written for the 1 _ state, the coefficients are given 
by C x . All the expressions obtained above have the same form, up to first-order in 
u) in the vicinity of u = 1, and up to the zeroth-order in 1/mg, as given by 

en = i - Q + <^ 2 >) - 1) , (54) 

A = lim E D = lim E B = C°, (55) 

rriQ-»oo mn—KXi 
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that is, 

e(l) = l, £'(!) = -i-Ii 2 <r 2 ). (56) 

The form for given by Eq. (|54p is derived from our semi-relativistic formulation 
and coincides with those derived from other considerations j23)J23,l22l we believe 
that our derivation is the most general and does not depend on any specific model. 
This form gives the lower bound —1/2 for £'(1) or the upper bound 1/2 for — £'(1), 
as shown in Refs. |24"|) . [27]) and |29|) . All the other form factors, £-(&>), £y(w), and 
£a ( w ) for i = 1 — 3, are similarly calculated and given by Eq. (|4ip up to this order, 
il)rn Q f. 

Using the values for the parameters obtained in Ref. fl"5"|) to calculate the mass 
spectra of heavy mesons, we can calculate £'(1) to first and second orders in 1/uiq 
as 

£'(1) = -1.44, (57) 

where we have used the values A = 0.752 GeV and (r 2 ^ = 5.009 GeV~ 2 , taken from 
Ref. I16|) . The values of £'(1) given in references are listed in Ref. 27). 

3.2. First order 

There is only one kind of the first-order correction, i.e., that from the wave func- 
tions. The first-order corrections to the form factors defined by Eq. f|39|) from the 
wave functions are straightforwardly calculated, although the calculation is cumber- 
some. Hence, here we do not repeat calculations similar to §3.11 and give only the 
final results below. 

We have calculated all four cases, 1-i) - 2-ii). The results show that there are no 
contributions from the first-order corrections of the wave functions to the form fac- 
tors. Although we attempted to obtain relativistically invariant results from knowl- 
edge of the rest frame by applying four different Lorentz-boost frames, the first-order 
corrections in all cases are not invariant, except one, the case 2-ii). The Breit frame 
with t' = x'° = y gives relativistic results that coincide with those of Ref. [21"]) . 
Hence, the following results could be a conjecture from our model, but we believe 
that our results are very plausible, since they agree with relativistic results which 
come from a semi-relativistic potential model. Brief derivations of these form factors 
are given in Appendix IE .21 

As stressed in the Introduction, our approach does not use fields, and hence 
there is no other kind of first order, i.e., the currents cannot be expanded in 1/rriQ 
in terms of the effective velocity-dependent fields, as in the HQET given in many 
papers as 

U = + 4 (£-7°K + ^7°7,?) K- 

Six form factors including zeroth-order terms for completeness with the decom- 
position of Neubert and RieckertPD 1 are given by 



1+ ( — + — ]px(w) 



(58a) 
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£-H = 




~ 2" 


CvM = 


1 + 




1 + 




1 


+ — 







2" 



1 



1 



1 



2w + l 



1 

1 1 



/ \ 1 
p 2 (uj) H 

1 / X 

H p 2 (w) 

rrir. 



P4H) 



1 



/54M 



£42 M 



1 



1 



cj + 1 m c 



[-A + {uj + l)p 3 (w) - P4(w)] f M 



1+ 2 



yl /u>- 1 1 



1 



, + — + — P2\U) 

uj + 1 m c mj, y m r . 



1 



P3(w) 



1 



W + 1 



p 4 (u) 



m 6 



where £(u) is given by Eq. ([33]) and 

p 3 (o;) = p 4 (^) 



-C^r 2 )^-!) 



0. 



(58b) 
(58c) 

(58d) 
(58e) 

(58f) 



(59a) 
(59b) 



Thus, there are 1/rriQ corrections to p\ and p 2 coming from phase factors of the wave 
functions, together with kinetic terms. (See Appendix IE. 21 ) 

§4. CKM matrix element 

4.1. B -» DID 

With the results given in § ^3. 1 1 and 13.21 we now calculate the CKM matrix 
element, \V cb \- We first evaluate the differential decay rate of the process B — > Dlv 
to extract the value \V C \\. Note that since we have neglected the isospin invariance 
in our formulation, there is no distinction between Bq and B~~ in the following. 
Similarly, £)W means either D or D*, which is equal to D° and D + or D*° and 
D* + , respectively, in our formulation. The expression for this differential decay rate 
is given by 



d£ _ _G 2 
~~ 48tt 3 

where we have defined 



3^ = T^T l^l 2 Mh (Mb + M D f (co 2 - l) 3 / 2 ^) 2 , 



1 + r 



M B 



(60) 



(61) 



Using the above expressions and the form factors obtained in the previous subsec- 
tions, we can evaluate values of the form factor Fd(u) at zero recoil and its first 
derivative. 
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i) Zeroth-order in 1/toq 

In this case, since = £y(^) = Ca 1 (^) = £,A 3 (&) = and £-(<-<->) = £a 2 ( w ) 

0, given by Eq. (f4l"]) . we have 



^b(i) = = i, ^b'(i) = r(i) = -i-44. 

ii) First-order in I/thq 

In this case, Eq. ([58]) . together with Eq. ([59]) gives 



(62) 



*b(l)=£+(l) 



1-r 
1 + r 



e_(i) = i.o7, ^b'(i) 



-0.875. 



(63) 



The combined recent experimental data of CLEO— ' and Belle^ lead to the 
following value of the product of the form factor and the CKM matrix element : 

Fd(1) \Vcb\ = 0.0414 ± 0.0064. 

Using our prediction Eq. (|63|). we find the value of CKM matrix element to be 

\V cb \ = 0.0387 ±0.0060. (64) 

4.2. B -> 

As noted by Neubertp^ the differential decay rate of the process B — > D*£u 
is the best quantity to extract the value |V^|. The expression for this differential 
decay rate is given by 

= 48^ |Fcfc|2 (MS " MD * )2 V^ 3 !^ + I)' 



1 + 



4uj 1 - 2wr* + r 
uj + 1 (1 - r*) 2 



*2 



(65) 



where we have defined 



*\2 



(1-r*) 



1 + 



4uj 1 — 2tur* + r 



+2 



uj + 1 {l-r*f 



2(l-2ujr* +r* 2 ) ( UA^Y + 



uj - 1 



■(w-lJ^M+rUM]} 



u; + 1 

2 



2(1 - 2wr* + r* z ) 1 + 



and 



UJ+1' 



iiiH 2 + [(a,-r*)-(o;-l)i? 2 (^ 



i2 2 (w) 



(66) 



(67) 



Similarly to §4. 1 1 one can evaluate values of the form factor Td*{^) at zero recoil 
and its first derivative. 
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Table III. Comparison of our model predictions for the values of the product Tr>*> (l)|Vci| with 
theoretical predictions and experimental dataP^ 

param. our Ref. ,33) 

Tp-\Vcb\ 0.0380(21) 0.0343(12) ~ 



params. CLECM^ 


BaBar 421 


Belled 


DELPHI 


, f 0.0376(3) (16) a 
>c * |Kc61 \ 0.0328(5)" 


f 0.0431(13)(18) a 
\ 0.0360(20) 6 


f 0.0354(19)(18) a 
\ 0.0315(12) b 


0.0377(ll)(19) a 



a fit using the form factor ( i a 1 parameterization of the model.^ 
b fit using the form factor predictions of the modelPS 



i) Zeroth-order in I/tuq 
In this case, we have 

^ D *(i) = = i, ^*'(i) = = -1.44. (68) 

ii) First-order in 1/mn 
In this case, we have 

T D * (1) = t(V) = 1, F D *'(l) = "1-09, (69a) 
Ri(l) = 1-45, R[(l) = -0.222, R 2 (l) = 0.942, i?^ 1 ) = 0.0286, (69b) 

where the first equation is consistent with the so-called Luke theoremp3 and we 
have used m c = 1.032 GeV, m b = 4.639 GeV, and r* = M D */M B = 0.3804 used in 
Ref. dZD and C 1 = 0.19022 in Eq. ([59a]) . The estimated values given by Eqs. ([63]) 
and ([69|) in our paper are consistent with other theoretical values, except for i?' 2 (l), 
whose values are given in Ref. [33]) . The value of i?2 (1) is one order of magnitude 
smaller than the other theoretical values, and there is yet no experimental data for 
it. 

The combined fit of the experimental data from CLEOpJ Belief and BaBar^ 
gives 

^b*(l) IV^I = 0.0380 ± 0.0021. 
Using our value of .Fd*(1) = 1, we have 

\V cb \ = 0.0380 ±0.0021. (70) 

§5. Summary and discussion 

In this paper, we have formulated a method for constructing Lorentz-invariant 
amplitudes by examining four kinds of Lorentz-boosted systems after normalizing 
their amplitudes. Although all four are not always consistent for higher orders in 
1/rriQ, they may be corrected by comparing with other methods, for instance, rela- 
tivistic and kinematic results. Despite of this defect, we believe that this formulation 
certainly provides a method to calculate reliable amplitudes from wave functions in 
the rest frame. 
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Table IV. Comparison of theoretical predictions of for the ratios -Ri(l) and -ffe(l) and their deriva- 
tives, R[(l) and R' 2 (l). 



Ref. 


flx(l) 


Ba(l) 


fli(l) 


TO 


our 


1.45 


-0.22 


0.94 


0.029 


El 


1.39 


-0.23 


0.92 


0.12 


IS} 


1.25 


-0.10 


0.81 


0.11 


Hz) 


1.27 


-0.12 


0.80 


0.11 


[20]) ,[48| 


1.35 


-0.22 


0.79 


0.15 


\M 


1.15 




0.94 




50) 


1.01(2) 




1.04(1) 





As an example of this formulation, we have calculated the semi-leptonic weak 
form factors up to first-order in 1/mg following a previous studypS 1 and using the 
results obtained in that paper. In this paper, we have found the following results. 

1. The Isgur-Wise function has the following form, up to first-order in I/thq and 
in (u — 1) : 

f (w) = 1 - ( i + \a 2 (r 2 ) ) (u - 1) , A = A u = lim E D = lim E B . 
Hence, we have 

ea) = i, e'(i) = -^-^ 2 <r 2 >. 

Here, since A depends only on the light quark mass and we treat only the heavy 
mesons D, D*, and B, which include only u and d light quarks with m u = ma, 
the subscript of A u expresses this fact. 

2. We find that there is no contribution to the six form factors from correction 
terms at first-order in thq for the rest-frame wave functions. That is, no terms, 
except for the first in Eq. (|27p . which include negative and positive energy states 
of a heavy quark contribute to the physical quantities. 

3. The first-order corrections are derived from phase factors of the wave functions 
and also given by kinetic terms, and there are no contributions from the first- 
order corrections to the wave functions. They are explicitly given by Eqs. (|58p 
and (|59p . That is, in the terminology of Neubert and Rieckert in Ref. |2T|) . we 
have 

fn(w) = p 2 (u) = ~C l A{r 2 ) (to - 1) , P3 (lo) = Pi {u) = 0. 

4. We have calculated the values for the form factor J-(uj) at zero recoil and/or 
their first derivatives up to first-order in 1/mg as 

^b(l) = 1.07, ^b'(l) = -0.875, 
JT D »(1)=£(1) = 1, F D J{\) = -1.09, 

fll(l) = 1.45, i?i(l) = -0.222, R 2 (l) = 0.942, R' 2 (l) = 0.0286. 

The first equations here were obtained by analyzing the B — ► Dlv process, 
and the second by analyzing B — > D*lD. These values are consistent with 
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experimental data as well as other theoretical estimates listed in Tables I and 
II of Ref. [33D- 

5. The above values can be used to estimate the CKM matrix element \V c b\, which 
we have obtained as 

\V cb \ = 0.0387 ± 0.0060 
for the B — ► Dlv process, and 

\V ch \ = 0.0380 ± 0.0021 

for the B — > D*£u process. These values are consistent with the value in 
PDGpJ 

\V cb \ = 0.0409 ± 0.0018 (exclusive). 

Here to obtain theoretical predictions for \V c f)\, we have ignored theoretical 
undertainties, although experimental errors are taken into account. 
We have developed a method to obtain relativistically invariant results using 
rest-frame wave functions, and to do this, four different Lorentz-boosted frames 
were adopted to check the validity of our results. The same form is obtained for the 
Isgur-Wise function in all four cases up to zeroth-order in 1/to.q and first-order in 
(oj — 1). However, the first-order corrections in 1/mn are not the same in all four 
cases. Only the case denoted by 2-ii) in the main text, i.e., in the Breit frame with 



JO 



t< 



given in Ref. [21]). 



y'° both for B and D^*\ gives results consistent with the relativistic ones 



Appendix A 

Schrddinger Equation 



The Schrodinger equation given by Eq. (J3j) is derived by calculating 
(X'; P x , {£'} \H\ X; P X {£}) = (X; P x {£} \V°\ X; P x {£}) , 



(A-l) 



with V° being the 0-th component of the four- momentum V^, and by varying with 
respect to ip x *- Here the Hamiltonian density is given by 



H-- 



d 3 x {q c] {x){a q ■ p g + (3 q m q )q c (x) + Q\x){a Q ■ p Q + f3 Q m Q )Q(x)} 



+ 



d 3 xd 3 y q c \x)f3 q O qi q%x)V ll {x-y^Q\y^l3 Q Qj Q{y). 

The lhs of Eq. (|ATj) can be approximated as 
(X';P X , {£'}\n\X;P x {£}) 



(A-2) 



X;P X {£} 



x I d 3 yq c J(t,x)q c a (t,x)Q}(t,y)Qp(t,y)H 



X;P X {£} 



~ / d 3 x / d 3 y(X';P x >{£'} Q}(t,y)q c a \t,x) 



x (0 



q c a (t, x)Qp(t, y)HQl(t, y)q^(t, x)qUt, x)Q s (t, y) X; P x {£} 
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~ Jd 3 xJ a?y^ x 7 a p((Q,x-y)-P x >) H a ^ s ^ s ((0,2 - y) ;P X ) 

x e - i ( p ^- p x')-y ; (A-3) 

where the number operators for q c and Q have been inserted, and the vacuum dom- 
inance has been used among the intermediate states. The rhs of Eq. (|A-1|) is given 
by 



(X'; P x , {£'} \V°\ X; P x {£}) = y/ M\ + P\ (X 1 ; P x , {£'} \ X; P x {£}) , (A-4) 
and thus we obtain the Schrodinger equation, Eq. ([3]). 

Appendix B 

Lorentz Boost and Normalization 



The derivations of Eqs. (|13p . (|15p and (|22p are similar to those given in 
Setting H = 1 in Eq. (|A-3p . we obtain 

(X;P X > {£'}\X;P X {£}) 

^Jd 3 xJ d 3 y ^ x * a/3 ((0,x-y);P x ,) tf Xaft ((0, x - g) ; P x ) e< p x- p x'* 

= (2vr) 3 5 3 (P x , - Px) J d 3 z tr ((0, z) ; P x ) tf x ((0, z) ; P x )] . (B-l) 
The lhs is given by 

(X;P X {£'} \X;P X {£}) = (2itf 2P Xo 5 3 (P x - P' x ) 5 e > e , (B-2) 

and hence 

J d 3 ztT[^ x \(0,z);P x )^ x ((0,z);P x )] ~2P XQ 5 ei (B-3) 

holds. 

We have to calculate the normalization of the wave function in the moving frame 
in two cases, i) t = x° = y° and ii) t' = x'° = y'°. To do so, we need to have a 
relation between the rest-frame (RF) and Lorentz-boosted (LB) wave functions, 
i) t = x° = y° and x'° ^ y'° (equal time in the rest frame) 
The relation between RF and LB is derived as follows. By definition, we have 



0\q c a (t,x)Qp(t,ff)\X- (M x ,0j {£} / 
= G~\Gp\ (0\q c J (x ,o ,x')Q s (y' o ,y")\X;Px {£}) . (B-4) 

Hence, we obtain 

^\ql{x' Q ^)Qp{y' Q ,y')\X-P x {£}) 

~ tf Xa p ((0, x' - g); P X ) exp {i[-M x t + (M x - m Q ) 7 2 /3(x 3 - y 3 )}) 
= G ai G pS (0 | q c ^t,x)Q 5 (t,y)\ X- (m x ,o) {£}) 

= G ai G ps i> Xl5 {x-y) e- iMx \ (B-5) 
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where use has been made of the approximation 

Q(3 (y' ^) * exp [-zm Q7 W° ~ x'°)] Qp (x' ^) . 
Rewriting Eq. (|B-5|) . we obtain the relation 

ric«f,(M;Px) * G ay Gp 8 ^s ((0,x±,7~V) ; (m x ,o)) e< M *- m 
Using this equation, we find 

d 3 x ((0, x);P x ) tf x a/3 ((0, x) ; P x ) 



(B-6) 



(B-7) 



7 / d\^*(x){G 2 ) ai {G 2 )^ XlS (x) 



7 / d x tr 



^(x)G^ x (x)G^ 



2MY 



(B-8) 



y /u and x° ^ y° (equal time in the moving frame) 



ii) t' = x 

Similarly to the case i), the relation between RF and LB is, in this case, given by 
;P X ) * G ay G ps ^ s ((0,f ± ,7x 3 ) ; (m x ,o)) e< M *~ m 

(B-9) 

Then, using Eq. (|B-9p we obtain 

d 3 x ^ ((0, 5) ; P X ) a/3 ((0, x) ; P x ) 

7- 1 / d?xi> x \p(x)(G 2 ) ai (G 2 )^ XlS {x) 



1 



1 / d 3 xtr 



^(x)G 2 i> x (x)G 



2T 



2M7. 



(B-10) 



Here use has been made of the following : 



1 

2 tr 



u_i(r) 
-i((T • n)f_i(r) 



t 



G 



2/ 



tr 





U-i(r) i(a ■ n)v-%(r) 

1 a 3 (3 



u-i(r) 
-i(a ■ n)v—i(r) 



1 a 3 P 

a 3 f3 1 



G 



2 T 



u_i(r) 
■ n)v-i(r) 



a 3 (3 1 



(B-ll) 



Appendix C 

Wave Functions 



The wave function is generally defined as 
(0\q c a (t,x) Qp(t,y)\X;P x {£}) 
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= (0 \q c a (0, V (x - y)) Qp(0, ((y-x))\ X; P x {£}) e~ iP ^ 

= A aP ms-y);Px)e- ip x-\ (c-i) 

with 

r» = (x^ + V y> M , ( + v = l. (C-2) 
Here, without loss of generality, we adopt 

C = o, 77 = 1. 

The Foldy-Wouthuysen-Tani transformation and charge conjugation are defined as 
t^FWT (?) = exp (w (p) 7q • pj = cos W + 7q • p sin W, (C-3a) 
tanW» = — *L_ E=y/p 2 +m* Q , (C-3b) 
^ = ^7q7q = -^ 1 - (C-3c) 

Here, p = pq is the initial momentum of the heavy quark, and C/fwt (p) operates 
on heavy quarks. Henceforth the color, iV c = 3, is ignored since the form of the 
wave function is the same for all colors. The quantity £ stands for a set of quantum 
numbers, j, m, and k, and 

*U*) = -( .-...^L * ) rfm(«), (C-4) 



r V -i v k{r) (n-a) ) 3 

where r = \x\, n = x/r, j is the total angular momentum of the meson, m is its z 
component, k is a quantum number which takes only values, k = ±j, ±(j + 1) and 7^ 

0, whose operator form is given by k = —f3 q i^J q • £ + 1^ . The scalar functions Uk{r) 

and Vk(r) are polynomials of the radial variable r and satisfy 

j dr(u 2 k (r) + vl(r))=l. (C-5) 

Here, y k m {fi) are functions of angles and spinors of the total angular momentum, 
j = I + s q + sq, with £ = —if x V, whose first few explicit forms are given by 

VoO = -?T=, Vim = -7T=° m 



/47T V47T 

Vim 



2 _ -1 3 



^n"-^!^ (C-6) 



These functions satisfy the follwing relation and normalization condition : 

Vjm = -{n-a) y k jm , (C-7) 



H/ 



dtf yj; m , V) 2/" m (^)) = 5 kk '5 jr 5 mm ,. (C-8) 
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The relative phases among y k j m given by Eq. (|C-6j) are determined so that they 
give the correct relative phases among the form factors, which are determined by 
Eqs. (|39aj) - (j39cp . From this point, whenever the trace of y^ m appears, it is under- 
stood to be in the sense of Eq. (|C-8p . 

The leading-order pseudoscalar state (0~) corresponds to £ = (k = — 1, j = m = 
0), and hence we have the wave function 

<*> - ( o iSi'M ) = o -i ) • < c ' 9 > 

with X = D or B. On the other hand, the leading-order vector state (1~) has the 
set of quantum numbers £ = (k = — 1, j = 1) and is given by 

i,' X0 (x) = s/M^ ( ) 

= M l -(° n »-><'■> ..W), (CIO) 
V 47r r \ • a)v-i(r) J 

with e m being a polarization vector (e 2 = —1) and X = D* or B*. 



Appendix D 

Matrix Elements 



The light anti-quark is treated as a spectator. That is, it does not interact with 
other particles, except for gluons represented by potentials in our model. The heavy 
quark has a weak vertex and its current is in general given by 



j(x,t) = Q\,Jt,x) O a pQ X p{t,x), 



(D-l) 



where O is a 4 x 4 matrix. By inserting the number operator of the anti-quark and 
by assuming vacuum dominance among intermediate states, we obtain the following 
formula for the matrix element of the above current, i.e., Eq. ()DTp . between heavy 
mesons : 



(X';P x ,,{£'}\j(t,x)\X;P x {t}) 



O a p(X' ] P x , {I 1 } 



Qx'afa*) / d 3 y q y(t,y)qC(t,y)Q X p(t,x) 



X;P X {£} 



j d z yQ a p (X';P X , {e'}\Q xla (t,x)q^(t,y) 
x (0\q^(t,y)Q X(i (t,x)\X;P x {£}) 

^(y-x;P xl ) V4 (y-x;Px) O 



or 



dry tr 



X';P X '{£'} j(0,0) X;P X {£} 



-i(P x -P x ,)-x 



(D-2) 



d 3 x tr 



^ x }{x;P x ,) rp x (x;P x ) O 



(D-3) 

where in Eq. ()D-2p we have assumed vacuum dominance among intermediate states 
in order to make some approximations. 
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Appendix E 

Semi-Leptonic Weak Form Factor 



In this appendix, we use the Foldy-Wouthuysen-Tani transformed wave func- 
tions, V'x FWT' instead of ip x , although the subscript FWT is omitted for simplicity. 

E.l. Zeroth order (Isgur-Wise function) 

Here we present the derivations of the Isgur-Wise function given in the main 
text in four cases. 
1-i) B rest-frame and t = x° = y° : 



2 ^ J d 3 x tr [Vdo ((0, x) ; P D ) ^ ((0, x) ; P) 
(l+u)V2m ^f d i xtT \^ Do{{ 0,x);P D )\ 2 

(l+uj) v / 2M B 2M D r 



xe 



-i(M D -m c )jVx 3 ,1 



2v^ 



2M(1 +w)7 3 / 2 



/ ^X \f D %aP (*) + - l)x 3 ^3^* 0Q/ 3 (*)} 



= I - <r 2 > - £LL| ^ x 3 (x) j ^ 0Q/3 (5) + 0(/? 4 ) 



(E-l) 



where we have used 



2M /" ,0 x 3 

= / cTx — tr 

4vr 7 2 



d 

dx 3 





u_i i{n-&)v-i 



_ 2M 
_ 2~M 

47T 

= -M, 



/ 



X' 

2 

x 3 3 



<i 3 x — tr 
d 3 x 



if 







r \ 





—i(n ■ a 




|(n 






a) 


r 


r 



4 dx 3 



tr 



(E-2) 



with M = ^M B M D . 

ae anc 

2v^ 



1-ii) 23 rest-frame and t' = x'° = y'° : 



(l+w)V2M B 2M D 7 



d 3 x g; 7 g; 5 Vdo 7 5 ((o,^,7^ 3 ) ; (m d ,o)) 



24 



T. Matsuki and K. Seo 



fd 3 x \ ^ D I ap (x ) + (7 - 1) x 3 ^3^ * a/3 (x) 



2M(l + a/) 7 1 /2 7 i™*^ ' w ' 8x 3 
xG a ,G^ e e 0l5 (x)e-^- m ^ x3 + 0(/3 4 ) 

= 1 - <r 2 > + ^-!=jd 3 xx 3 (Jgrilap (*)) (*0 + O(^) 

= 1 - i/? 2 u,^ <r 2 > - ^ + 0(/3 4 ). (E-3) 

2-i) Breit frame and t = x° = y° : 

/ d 3 x tr U£ f ((0, x) ; Pjj) Vfjo ((0. I p b) 



/f^ 3 * trl^oaO.fJjPr,)! 2 Jd 3 x tr |^ ((0, £) ; P fl )| 



d 3 x 



V / 2M B7 3 2M D7 3 



2M 7 3 



2 



= 7~ 2 - I (f) 7 2 (^ + ^d) 2 <r 2 > + 0(/? 4 ). (E-4) 
2-ii) Breit frame and t' = x'° = y'° : 

^ = jmkm=, I d3x <vw£5* ; (m.,o)) 



2M 7 

= 7" 2 " ^ (f ) 2 7' 2 (^ + Pd) 2 <r 2 > + 0(/3 4 ). (E-5) 
E.2. Firsi order 

In this appendix, we show how to calculate the first-order corrections to the 
form factors by using Eqs. (f25j) - (f2"7"|) and (f29j) and the equations in Appendix O 
Only the case 2-ii) gives six form factors that are consistent with those of Neubert 
and RieckerlP among the four different Lorentz frames considered in the previous 
subsection : 
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2M 7 
1 

1 



l_ f 

f 7 2 J ' 



i4 t 



1 



5 + ^ 2 <r 2 >]U'-L) 



1 ,, 1 ^ ^ 
— 7 • Pc 7 • Pfc 



2 



1 ^ 1 



+ 



l--IC 1 <r 2 >(u,-l)(^- + ^- 

l-^<r 2 >(.,-l)(^- + ^- 
6 \m c nib 



1 - 



( W -l) 



(E-6) 



where for X = B and/or D are defined by Eqs. (fIBj) and (|C-3l) and M 
y/M B M D . Here, yl is denned by 



A 



Ep + E B = (Ed - m c ) + (E B - m b ) 
2 2 



(E-7) 



which comes from the expansion of the phase factor exp (-2iA^z 3 ). Here, we have 



C° = A = A u and C 1 is defined in Eg. ([53]) and depends only on 
Comparing the result with Eq. (|58aj) . we obtain 

1 1 



m n = m,, 



^) = l-[- + -A* 



( W -l) 



m d . 



(E-8) 



!+(— + — )pi(wj 



£M, Pl (u;) = --iC 1 <r 2 >( W -l). (E-9) 



Similarly, the five other form factors and Pi(to) are given by the following : 
1 



2M 7 2 § 



<i 3 z tr ^^^^^t/^^-^^G*^) e" 



-2iA§z 3 



1 



2M 7 2 



d 3 *tr(^ ^U^(a 3 + ^ 



1 



1 



1 3 
Pc + Pb 



A 

"2 +Pi( u} ) 



m h J 



cm, 



(E-10) 



(E-ll) 



frM = -= 



1 



2M*7 2 I(-^ 2 *) 



2M* 7 f(-ie 2 *) 



d 3 z tr ^^(^^(^^b 10 " 1 "^ *^*) e" 
d 3 z tr (^VWioC*) 
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1 3 1 3 
—Pt Pb 



7 °E 2 



1 + 



A' ( 1 



+ 



1 



2 \m c rrif, 



yl / 1 1 

1 + IT — + — I 

2 \m c mi, 

p 2 (v) = pi(u), p 4 (cj) = 0, 



— + — ) pi(w) 



CM, 



(E-12) 
(E-13) 



where Z 1 



a 
Here, we define 



a 



, e is a static polarization vector for D*, and M = \JMbMd* 
, _ E D * + Eg _ (E D * - m c ) + (E B - m b ) 



C 



c 1 



= C ° + 2^ + 2^ + °( 1 /^ 2 )' 
with C ' = C° = /l = A u , 



(E-14) 



2M 7 2 
1 



2M Ve 1 * 



d 3 z tr 



(^d* V)V4o(*) 



1 3 1 3 
— 



7°^ 



-2iA'S-z 3 



1 + 



A'u-l ( 1 



+ 



1 



2 cj + 1 \ m c wi& 



1 + 



ylw - 1 / 1 



+ — ] P2(w) 



2 W + 1 \ m c ^6 / V m c Wife 

p 2 (w) = piH, p 4 (w) = o, 



en, 



(E-15) 
(E-16) 



~Ui H + 6t 2 (w) + £a 3 M 

= - 2 m* 7 ^ £ 3, tr (^-oW*,,^ (l + 5 



1 3 1 3 
— Pi + — P6 



e -2i/l'fz 3 



yl' / 1 



M 

2 ^m c m ft y 



(E-17) 
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1 f d 3 ztr^ D , \z)^ Bo (z)^ 1 



Pc Pb 

m c m b 



757° ) ) e" 2 ^ 3 



1 " ^) fa)- (E-18) 



w + 1, 
From these, we derive 



CiaH = fa), Pz{oj) = p A {ui) = 0, (E49) 



u + 1 m c 

1 1 \ / 1 1 , 

2 \ w + 1 m c mj / V m c m& 



(E-20) 



p 2 (w) = Pi(^), P3(^) = AtM = 0. (E-21) 

Therefore, all the expressions for pi(u>) included in £j(a>) are consistent up to first 
order in 1/itiq. Even though the first-order corrections to the wave functions are 
included as given by Eq. (|27p. their contributions vanish, because of the matrix 
structure after taking the trace over indices. From the above equations, we can easily 
reproduce our results given by Eq. ([58]) . together with Eq. ([59|) . 
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